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Randomness-Driven Quantum Phase Transition in Bond- Alternating Haldane Chain 
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The effect of bond randomness on the spin-gapped ground state of the spin-f bond-alternating 
antiferromagnetic Heisenberg chain is discussed. By using the loop cluster quantum Monte 
Carlo method, we investigate the stability of topological order in terms of the recently proposed 
twist order parameter [M. Nakamura and S. Todo: Phys. Rev. Lett. 89 (2002) 077204]. It is 
observed that the dimer phases as well as the Haldane phase of the spin-1 Heisenberg chain 
are robust against a weak randomness, though the valence-bond-solid-like topological order in 
the latter phase is destroyed by introducing a disorder stronger than the critical value. 

KEYWORDS: Haldane chain, bond randomness, quantum phase transitions, random-singlet phase, quan- 
tum Monte Carlo, loop algorithm, topological order, twist order parameter 
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Disorder effects on low-dimensional quantum magnets 
have been investigated extensively in recent theoretical 
studies. In particular, impurity effects on spin-gapped 
Heisenberg antiferromagnets 1 have aroused much inter- 
est in relation to the impurity-induced antiferromagnetic 
(AF) long-range order observed experimentally in real 
materials. 2 It has been established by recent numerical 
simulations 3,4 that in two dimensions or higher, there 
are two classes of disorder, that affect spin-gapped states 
in essentially different ways. Site dilution and bond di- 
lution are representatives of each class. The former in- 
duces localized moments around impurity sites. There 
exist strong correlations between such effective spins re- 
taining the staggeredness with respect to the original lat- 
tice, and therefore the AF long-range order emerges by 
an infinitesimal concentration of dilution. In the bond- 
dilution case, on the other hand, localized moments are 
always induced in pairs and they form a singlet again by 
AF interactions through the two- or three-dimensional 
shortest paths as long as the concentration of bond dilu- 
tion is smaller than a finite critical value. 

In one dimension, since quantum fluctuations are much 
stronger than those in higher-dimensional systems, novel 
quantum critical phenomena are observed under disor- 
der at the magnitude of coupling constants (bond ran- 
domness) . Theoretically, the decimation renormalization 
group (DRG) approaches have achieved great success in 
predicting rich physics, such as the random-singlet (RS) 
phase for spin-i chains. 5-7 Recently, this technique has 
been extended to higher-spin cases, 8-11 where two of the 
main debates are on the robustness of the Haldane gap 12 
against disorder and on the presence of the spin-1 RS 
phase. A number of numerical studies have also been 
carried out 13-16 to establish a quantitative phase dia- 
gram. However, this problem has not been sufficiently 
clarified yet. One of the main difficulties in simulating 
random quantum systems is the extremely wide energy 
scale that has to be taken into account. Another difficulty 
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is the lack of appropriate physical quantity for effectively 
discussing randomness-driven critical behavior. 

In this Letter, we report the results of our quan- 
tum Monte Carlo (QMC) simulation on the bond- 
alternating Haldane chain with bond randomness. By 
using the recently proposed twist order parameter 17 to- 
gether with a novel numerical technique for simulating 
the ground state in the framework of the loop cluster 
QMC method, 18-20 we show that the difficulties men- 
tioned above can be overcome. Thus, we successfully es- 
tablish the quantitative ground-state phase diagram. 

We start with the following Hamiltonian for the AF 
Heisenberg chain 



(1) 



Here, Sj is the spin-1 operator at site j and L the system 
size; periodic boundary conditions are imposed. 

For the bond-alternating model without disorder, 
where the coupling constants {Jj} are given by Jj = 
1 — (— 1) J 6 parameterized by the strength of bond alter- 
nation 8, its ground state has been discussed in terms 
of the valence-bond solid (VBS) picture. 21 For the spin 
size S, the pattern of the valence bonds (to, n), where 
m (n = 2S — to) denotes the number of effective sin- 
glet pairs on odd (even) bonds, changes from (0, 2S) to 
(2S, 0) successively as 8 is increased from —1 to 1, indi- 
cating the existence of 2S quantum phase transitions. 22 
Each VBS state has a topological hidden order, which is 
characterized by the string order parameter. 23 

On the other hand, Affleck and Lieb studied Haldane's 
conjecture on the basis of the Lieb-Schultz-Mattis (LSM) 
argument. 24 Although the association between the VBS 
picture and the LSM argument has not been fully under- 
stood for a long time, Nakamura and Todo have recently 
shown that the ground-state expected value of the uni- 
tary operator appearing in the LSM argument, given by 



z L = (exp[i^ J2jS?]) , 

3=1 



(2) 
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functions as an order parameter, which characterizes the 
VBS states. 17 The unitary operator in eq. (2) rotates the 
spins about the z-axis with the relative rotation angle 
2n/L; thus, it generates a low-lying excited state with 
an excitation energy of C(L _1 ). Since the twist order 
parameter (2) measures the overlap between the ground 
state and such a twisted excited state, \zl\ ^ 1 in the 
thermodynamic limit evidences the existence of gaplcss 
low-lying excitations or a degeneracy in the ground state. 
Furthermore, it is shown that in the (m, n) VBS phase, 
Zl converges to (— l) m for L — > oo. We will see below 
that the twist order parameter works fairly well even in 
the presence of disorder. 

In what follows, we consider two different random dis- 
tributions for the couplings {Jj} in eq. (1). The first one 
is the uniform distribution, where the coupling constants 
are distributed uniformly according to 



\l/2W if \Jj - 1 + {-l) j 8\ < W 



otherwise. 



(3) 



Here < W < 1 — \6\ must be fulfilled, otherwise ferro- 
magnetic bonds could appear in the system. The second 
distribution is given by 



(4) 



jj — l± y ±j uj uj 

with the quenched random numbers tj obeying the 
power-law distribution 5,6 ' 25 



'R-Hj 1+1/R if < *j < 1 







otherwise, 



(5) 



with a non-negative parameter R, where the R — > limit 
corresponds to the nonrandom case {tj = 1 for all j). 
Note that at 5 = the uniform distribution (3) with W — 
1 is equivalent to the power-law distribution [eqs. (4) and 
(5)] with R — 1 besides a trivial scaling factor; Jj's are 
distributed uniformly between and a finite cutoff. In 
the following simulations, we take the random average 
over 1000 samples for each parameter set. 

The present model (1) can be simulated efficiently by 
the loop cluster QMC method 18 ' 19 even in the presence 
of randomness. However, it should be pointed out that 
the loop cluster method, which is based on the Suzuki- 
Trotter path-integral representation, works indeed at a 
finite temperature. Since the ground-state properties are 
only our main concern in the present study, an effective 
extrapolation scheme, which we will explain below, for 
taking the zero-temperature limit is essential. 

We notice the fact that the ground state of the nearest- 
neighbor AF Heisenberg chain of finite and even number 
of spins is singlet, and there is a finite gap above the 
ground state. In the path-integral representation, the in- 
verse of the gap is given by the correlation length along 
the imaginary-time axis. Since loop size is directly re- 
lated to the correlation length in real space as well as in 
the imaginary-time direction, 18 ' 19 the system cannot dis- 
tinguish whether the temperature is finite or zero, if no 
loops wrap around the lattice in the imaginary-time di- 
rection. In other words, the winding number of the loops 
in the imaginary-time direction can be used as a good 
measure of the convergence to the ground state. 
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Fig. 1. (^-dependence of twist order parameter z^ for spin-i sys- 
tem with R = 0.5. At S = 0, Zl is zero irrespective of system 
size, while it converges to ±1 for 8 ^ 0. 



Although there exist several means of implementing 
the above idea as a ground-state QMC algorithm, 20 ' 26 
we employ the following in the present study. We start 
with a certain temperature. During the thermalization 
Monte Carlo sweeps, the winding number of the loops is 
monitored. If one or more loops wrap around the system 
in the imaginary-time direction, we double the inverse 
temperature. This procedure will automatically adjust 
the simulation temperature so that the system will be at 
the ground state effectively. 

Before jumping to the spin-1 system, we discuss briefly 
the phase diagram of the spin-i system, for which the 
effects of disorder on this system have been well es- 
tablished. The ground state of the non-bond-alternating 
spin-i chain without disorder is critical. By introducing 
infinitesimal randomness, the system is driven to the RS 
phase, where there is also no excitation gap, but the cor- 
relation function decays with an exponent different from 
that of the nonrandom system. 6 The RS phase is charac- 
terized by an infinite dynamical exponent, i.e., a logarith- 
mic scaling of the length and energy scales. As a result, 
the uniform susceptibility diverges as \ ~ l/7Tog 2 T at 
low temperatures. 27 

The RS phase is unstable against bond alternation. 
The real-space correlation becomes short-ranged imme- 
diately, though the spin gap remains absent up to a finite 
strength of bond alternation. 7 ' 27 This phase is referred 
to as the quantum Griffiths (QG) phase, where the uni- 
form susceptibility obeys the power law {\ ~ T~ 7 ) at 
low temperatures with a nonuniversal exponent 7 vary- 
ing with 5. 

In Fig. 1, the twist order parameter is plotted as a 
function of 5 for the spin-i chain with R = 0.5 (power- 
law distribution). The twist order parameters with dif- 
ferent system sizes clearly cross at 8 = 0. Note that in 
the random system, the translational and parity symme- 
tries are both broken in each sample, and thus zt, does 
not necessarily become zero at 5 = 0. However, one sees 
in Fig. 1 that the symmetries are restored after the ran- 
dom average is taken. For a nonzero 5, the twist order 
parameter rapidly converges to ±1, though gapless QG 
phases extend on the both sides of the RS point. 7,27 The 
present results demonstrate clearly that the twist order 
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Fig. 2. Ground-state phase diagram of spin-1 chain with uniform 
random-bond distribution. Along the <5 = line, the Haldanc 
phase survives up to W = 1. 



Fig. 3. Scaling plot of distribution function of local susceptibility 
at RS point (R, S) = (0.5, 0.2). Inset: 5-depcndence of twist order 
parameter for R = 0.5 (crosses, squares, and triangles for L = 8, 
16 and 32, respectively.) 



parameter zl is not affected by QG singularity, and thus 
it is an effective tool for analyzing RS criticality. 

In contrast to the spin-i chain, the non-bond- 
alternating spin-1 Haldane system without disorder has 
a finite gap and a finite correlation length. 12 In the pre- 
vious DRG studies, 8,9 it is predicted that the Haldane 
state is stable against a weak disorder, while there occurs 
a quantum phase transition to the spin-1 RS phase at 
a critical strength of randomness. In the previous QMC 
analysis 15 of the model with a uniform random-bond dis- 
tribution (3), in which the uniform susceptibility and the 
string order parameter were mainly investigated along 
the 5 — line, it was concluded that a quantum phase 
transition occurs at W ~ 0.95 from the Haldane phase 
to the RS phase. In the present calculation, however, the 
twist order parameter decreases with increasing system 
size in the entire range of W (0 < W < 1), and tends 
to converge to -1 without showing any crossing, which 
indicates that the Haldane [(1,1) VBS] phase is stable in 
the entire range of W. 

This can be seen more clearly in the S-W phase di- 
agram shown in Fig. 2. The phase boundaries are ob- 
tained from the crossing point of the twist order param- 
eter with different system sizes (L = 8 • • • 64) . For small 
(5's, where the Haldane phase existing at \S\ < 0.25997(3) 
for W — 17 decreases gradually, the phase diagram 
[Fig. 2] agrees qualitatively with the one predicted by the 
DRG analysis. 10 However, the phase boundary between 
the Haldane (1,1) and the dimer (2,0) phases (solid line) 
merges with the parameter boundary S + W = 1 (dashed 
line) at S ~ 0.1, and does not reach 5 = even at W = 1, 
indicating that there is no spin-1 RS phase in the model 
with the uniform random-bond distribution. 

Next, we examine the other random-bond distribution, 
i.e., the power-law distribution [eqs. (4) and (5)]. As al- 
ready mentioned, the power-law distribution with R = 1 
is equivalent to the uniform one with W = 1; thus, it is 
expected that the Haldane phase is stable at least up to 
R = 1 also for the former case. However, for the power- 
law distribution, one can consider a further strong disor- 
der (R > 1), i.e., a wider distribution on the logarithmic 
scale, by which the Haldane phase might be broken. 25 

In the inset of Fig. 3, the twist order parameter is plot- 



ted as a function of S in the weak randomness regime 
(R = 0.5). For S > 0.2, the twist order parameter in- 
creases as the system size increases and tends to con- 
verge to +1. We identify this phase as the dimer (2,0) 
phase. On the other hand, zl tends to converge to -1 for 
S < 0.2, indicating the Haldanc (1,1) phase. 

At the crossing point 5 = 0.20(1), a quantum phase 
transition occurs, and the transition is expected to be- 
long to the spin-i RS universality class. 10 To confirm 
this prediction, we measured the distribution of the local 
susceptibility 

Xioc,,-=/?K 2 >= [ dT(S*(O)S*(r)) (6) 
Jo 

at the critical point (R, 8) = (0.5, 0.2). As seen in Fig. 3, 
the distribution function of the logarithm of the local 
susceptibility is scaled fairly well by assuming a loga- 
rithmic scaling form, P(logxioc) — /(logXioc/-^) with 
tp = 0.42. This is consistent with the previous DRG 
prediction for the RS phase, 6 though the exponent ip 
is slightly smaller than the predicted value (ip = 1/2). 
This is additional support to the validity of applying 
the twist order parameter to randomness-driven quan- 
tum phase transitions. Repeating similar analyses, we 
obtain the entire S-R phase diagram of the random Hal- 
dane chain with the power-law distribution [Fig. 4]. 

Although the phase diagram for the power-law distri- 
bution is similar to Fig. 2 for small R's, the overall shape 
of the phase boundaries indicates the existence of a mul- 
ticritical point, where two critical lines merge with each 
other at a finite R. To locate the multicritical point, we 
calculate zl for several system sizes (L = 16 • • • 64) along 
the 6 = line. The results for 0.9 < R < 1.2 is shown in 
Fig. 5, where the data with different system sizes clearly 
cross at R c ~ 1.05. Thus, we conclude that there exists 
a multicritical point at (R,S) = (1.05,0), which is indi- 
cated by the solid square in Fig. 4. Below the multicrit- 
ical point, the Haldane phase survives, though the spin 
gap vanishes at a certain R (< R c ), where a crossover 
from the gapped Haldane phase to the gapless Haldane 
(or QG) phase occurs. In the case of a uniform distribu- 
tion (3), the crossover is observed at W ~ 0.7, 27 though 
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Fig. 4. Ground-state phase diagram of spin-1 chain with power- 
law random-bond distribution. The multicritical point is indi- 
cated by the filled square. 
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Fig. 5. -R-dependence of twist order parameter zl of spin-1 sys- 
tem with 8 = 0. The data with different system sizes cross at 
R ~ 1.05. 



we have not yet examined it for the power-law distribu- 
tion. For R > i? c , on the other hand, the twist order 
parameter is expected to converge to a nontrivial value 
(i.e., Zl ^ ±1) in the thermodynamic limit, where the 
spin-1 RS phase is realized. 8-10 

To summarize, we reported the results of our QMC 
simulations on the bond-alternating random Haldane 
chain. By introducing the ground-state loop cluster QMC 
method and the twist order parameter, we have success- 
fully calculated the ground-state phase diagram. In par- 
ticular, we demonstrated that the twist order parameter, 
introduced originally for pure spin chains, is also effec- 
tive for random spin chains. Indeed, it is shown that the 
behavior of the twist order parameter observed in the 
present study can be discussed more directly in terms of 
the numerical DRG approach, in which one can calcu- 
late the topological order parameter for an approximate 
VBS-like ground state explicitly. 28 

For the uniform distribution, the present result, i.e., 
the absence of the spin-1 RS phase, does not agree with 
the previous finite-temperature QMC result, in which 
a multicritical point was suggested. 15 The possible rea- 
son for this disagreement is that the finite-temperature 
QMC method might easily fail to take into account rare 



and low-energy-scale but very strong correlations, which 
are essential in random spin systems. In contrast, in the 
present ground-state algorithm, simulation temperature 
is automatically adjusted according to the magnitude 
of the gap of each random sample, so that the physi- 
cal quantity at the zero temperature is calculated at an 
optimal cost. This algorithm is useful in simulating not 
only random systems but also those without disorder. 20 

For the power-law distribution, on the other hand, we 
established a phase diagram with a multicritical point, 
whose location was also determined accurately using the 
twist order parameter. The present phase diagram agrees 
qualitatively with the recent DRG prediction, 10 though 
the numerical confirmation of spin-1 RS criticality, which 
is expected to realize in the strong disorder regime (R > 
i? c ), still remains as a future problem. 
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